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Abstract 

We present a new model for credit index derivatives, in the top-down approach. This model 
has a dynamic loss intensity process with volatility and jumps and can include counterparty risk. 
It handles CDS, CDO tranches, N th -to-default and index swaptions. Using properties of affinc 
models, we derive closed formulas for the pricing of index CDS, CDO tranches and N th -to-default. 
For index swaptions, we give an exact pricing and an approximate faster method. We finally show 
calibration results on 2009 market data. 
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1 Introduction 



Credit derivatives allow to transfer some credit risk: an investor can get or sell protection on a debt 
issuer or on a basket of different names. Most popular credit basket derivatives are Credit Default 
Swaps (CDS), where one gets protection on a whole basket against some premium paid on a regular 
basis. For CDO tranches, the protection does not cover losses below some trigger or above a second 
one. In N th -to-default, the protection concerns precisely the N th default on the basket. An other kind 
of derivatives are CDS options (swaptions): one gets the right to enter at a future date in a CDS at a 
given premium rate. Depending on the market price of protection at the option maturity, the option 
is exercised or not. 

There are different risks in credit derivatives. The risk of default should be the more obvious. A 
second kind of risk are movements in the market price of CDS. When the derivative is written on a 
basket, risks comes from default risk and spread risk on under lyings and also depends on correlation 
between names. 

There is a first distinction between credit models. Some models are static: the intensity of default 
of a name has no stochastic behaviour with respect to time evolution. It is sufficient to price CDS 
or CDO tranches, but such models cannot handle options: the spread risk is not taken into account. 
Dynamic models incorporate this risk, generally by making the default intensity stochastic. 

There is an other splitting between credit basket models. In bottom-up models, credit basket 
derivatives prices and risks are reconstructed from individual constituents and their correlations 
|LiOO| IDC01] . An alternative approach describes directly the risks on the basket |EGG09| [Sch06, 
IBPT071 IAH07| . It is sufficient if one wants to hedge derivatives on a credit index with the index itself 
or other derivatives and not with individual components. If one is concerned with individual names 
risk, a second step is necessary, such as the random thinning [GGD09 . 

We introduce a model which belongs to the second category: we model directly defaults on the 
basket. Our model is also dynamic: the intensity of default is a dynamic, stochastic variable. More 
precisely, the intensity of default of the basket is modeled by a Cox-Ingersoll-Ross process with jumps. 
A diffusion term describes short movements of the credit index spread. Jumps intend to take into 
account possible crises where the perceived credit quality of the basket decreases by a large amount. 
After crises, the intensity can relax to a lower value. We also introduce a probability for the default 
of the whole basket. Our model can handle counterparty risk: the counterparty can be correlated to 
the basket or can even be one of the basket names. 

We define first a model on an infinite pool of names where default events occur according to a 
stochastic default intensity, without any bound on the total number of defaults. In a second time we 
use this infinite pool model to model defaults on a finite basket. Using properties of affine models 
[DPSOO^, we get closed formulas for the (conditional) expected loss or tranches losses. Usual credit 
derivatives can be priced exactly: CDS, CDO tranches, N th -to-default, index CDS options. It takes 
reasonable time for most instruments. However for swaptions the exact pricing can be slower so we 
introduce some approximations which reduces highly the computation time. 

We describe our model in section [2] and the pricing of usual credit derivatives in section [3] We 
finally show numerical results of calibration on market data in section |4j 
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2 Model 



2.1 Infinite pool 

2.1.1 Description of the model 

The credit portfolio in modeled in a top-down approach: the portfolio loss is considered directly. 
Probabilities are considered to be risk-neutral probabilities. We consider first a portfolio of infinite 
size. Its loss Lt is a stochastic process with jumps I dJ which occur with an intensity At and size i. 

dL t = I dJ 

where dJ are jumps of size 1. The recovery rates R of defaults are independent and identically 
distributed. We denote by dvj(l) the distribution of the jump size I = 1 — R. In parallel, we have a 
stochastic process Nt counting the number of default. It is similar to Lt but with zero recovery rate: 

dN t = dJ . 

We want to build a model for the default intensity process with the following features: 

• a diffusion process describing small fluctuations of the CDS index spread; 

• positive jumps of the intensity which can produce episodes of crisis with many defaults on some 
timescale; 

• a mean reversion such that crises can relax. 

Mathematically, the intensity At will follow a stochastic process like 

d\ t = k(Aoo - X t )dt + a{X t )dW + dX 

where dW is a Brownian motion and dX a jump process with intensity 7. 

It appears that this model has analytic solutions if we make some further assumptions on the 
diffusion and jump process: 

• the volatility is proportional to the square root of A: cr(\t) = o~\fXt, as in a Cox-Ingersoll-Ross 
process; 

• the jump size are independent and identically distributed along a Gamma distribution of pa- 
rameters n and 9 where n is an integer which describes the shape of the distribution and the 
mean size of jumps is (n + 1)0: 

x n e -x/e 

dv x{x) = nWn+1 l x>0 ; (1) 

• all processes are independent of each other (except the fact that At is the intensity of jumps 
dJ). 

The justification of this set of assumption is computational tractability. This choice encompasses the 
features we want to model without restricting the model too much. In fact one feature we might want 
to include is missing: here jumps in the hazard rate are purely exogenous, there is no default contagion, 
as in [EGG09 . Including a jump term proportional to dLt for the intensity process requires numerical 
integration of differential equations, which is tractable but much slower than analytical integration. A 
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few tests do not indicate that including default contagion changes the loss distributions and therefore 
the CDO prices. 

The default intensity process At is thus supposed to follow a Cox-Ingersoll-Ross process with 
additional jumps with intensity 7 and size distributed along a Gamma distribution (equation (U])): 



dAt = k(Aoo - X t )dt + a^X t dW + dX . 

Parameters k, Aqo, a, n, 7 and 9 can be taken to be constant or piecewise constant. The model 
parametrization contains in addition the initial default intensity Ao- 

Experience shows that calibration of the model without a probability of default of all names tends 
to incorporate it by setting a very high jump size in dX. We therefore introduce explicitly such a 
possibility. We introduce a variable Q which jumps of one unit when all names default together, with 
intensity aXt + 0. Q 7^ models a default of the whole portfolio. Parameters a and (3 control the 
probability of such an event and its correlation to the short-term credit quality of the basket modeled 
by A. a gives a contribution independent of the current state of the system whereas (3 gives a higher 
probability of total default during crises. 

In order to handle counterparty risk, we introduce finally a variable R which jumps of one unit 
when the counterparty defaults. It jumps with probability £ when there is a jump in dJ and addi- 
tionally with intensity QXt + t] at any time. When the counterparty belongs to the set of index names, 
£ and r\ are set to zero and £ controls the probability that a default concerns the counterparty. When 
the counterparty is not one of the names, £ is set to zero, £ and r\ control the intensity of a default 
of the counterparty and its correlation with defaults in the basket. 

The intensity process jump distribution dvx(x) given in ([TJ can be more general than a pure 
gamma distribution by considering the sum of several gamma distributions with different parameters 
7®, re® and 9^. If 9 is kept constant, one can for example take a distribution of the form P(x)e~ x ^ e 
where P is any polynomial. For clarity we write formulas for only one gamma distribution; the 
generalization is straightforward. 

2.1.2 Characteristic function 
2.1.2.1 Theorem 

Theorem 1. For constant parameters, the characteristic function at date t in the four variables Ly, 
Nt, At and Qt for complex variables u, v, w and x is 

E| ^ u Lt+vN t +w\ t +xQ t +vRt\ _ e A(t,T;u,v,w,y)+B(t,T;u,v,w,y)Xt+uL t +vN t +xQt+yR t 

witl^\ 

B+ - B_ 
B = B_ + ^-^- 

x 

k ± v 7 ^ 2 + 2a 2 ip 



B± 



= W ~ B- c _ y/K 2 +2a 2^ (T _ t) 

w — B + 

tfj = l- e v <t>j(u) [1 - £(1 - e y )] + a(l - e x ) + £(1 - e y ) 



To simplify notation we drop the explicit dependency of functions A, B, B±, x an d tp in variables t, T, u, v, w, x 
and y. 
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A 



\o\B+ - B.) 



k\ 00 B± + 7 
1 

(1 - 6B±) n + l Ow 



(l-6B±) n + l 



In 



B-B± 
w-B ±i 



Al l_ 

1 ^k(l-6B+ 



k (i — pr>-|- 

-^-^(r-o-^i-e^tr-t) . 



in- fc+1 



1 



(l-0B) fc (1-0™)* 



4>j(z) is the characteristic function of the individual loss size: 

<j>j{z) = [ e zl dvj{l) . 



The proof of this theorem and its generalization to piecewise constant parameters is described in 
the following. 

2.1.2.2 Differential equations 

Proposition 1. The model we have described is an affine model in the three variables Lt, Nt and Xt: 
its characteristic function can be written in the form 



jg / ^uLt+vNt+wXt+xQt+vRi 



^A(t,T;u,v,w,x,y)+B(t,T\u,v,w,x,y)\t+C(t,T;u,v,w,x,y)Lt+D(t,T;u,v,w,x,y)Nt+E(t,T\u,v,w,x,y)Qt+F(t,T\u,v,w,x,y)Rt 



(2) 



where A, B, C , D, E and F are complex functions of the real parameters t and T and the complex 
parameters u, v, w, x and y which satisfy differential equations and boundary conditions 

pb -dis- 



proof. Let suppose that equation @ holds. The variation of E t l^L T +vN T +w\ T +xQ T + y R T ^ when t 
goes to t + dt can be written as 



dE [ e uL T+vN T +w\T+xQ T +yRT 



d t Adt + d t Bdt X t + (e 



d t Edt Q t + d t FdtR t + ( e C<il t +DdN t +FdQ t _ A + / e Ed Qt _ 1 



1 +d t CdtL t + d t DdtN t 



D A+B\t+CL t +DNt+EQt+FRt 



(3) 



where 



1 = BK(\ O0 -X t )dt + Bay / X t dW+-B 2 a 2 X t dt+ [e BdX -1 



CdL t +DdN t +FdRt _ j = e (Cj+D)d,J+FdR _ 1 



On the other hand, we know from the trivial equation 



Et+dt ^e u LT+vN T +w\ T +xQ T +yR T '^ _ jg^ ^uLT+vNT+wXT+xQT+yRrj 



that 



E* 



d E ( e u LT+vN T +w\ T +xQ T +yRT 
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Therefore the expected value of (pi) must vanish: 



+ t B it A t + B .(A M - A«)d t + 'bV^ + Td t / (e- - 1) *»(.) 
+ tCit Z, + *D* JV, + e tEit Qt + 6,F it Rt + A, dt (e° / e^ ) + (1 - 0] ~ l) 



(aA t + /3)dt (e E - l) + (CA t + r/)dt (e F - l) 



Forgetting the last exponential, this equation is affine in Lf, Nt, Xt, Qt and Rt- As this must be 
true for any value of stochastic variables L t , N t , Xt, Qt and Rt, this equation holds if and only if the 
following equations are satisfied 



d t A(t,T;u,v,w,x,y) -- 

d t B{t,T;u,v,w,x,y) = 

d t C(t,T;u,v,w,x,y) = 

d t D(t,T;u,v,w,x,y) = 

d t E(t,T;u,v,w,x,y) = 

d t F(t,T;u,v,w,x,y) = 

with the boundary conditions at t 



-kX^B - 7 (ct> x (B) - 1) - (e E - 1) - 77 (e F - l) 

«B " \* 2 B 2 - (e D MC) & F + (1 - 0] - 1) 

-a (e E -l)-C(e F -l) 







. 

T 



A(T,T,u,v,w,y) = 

B(T,T,u,v,w,y) = w 

C(T,T,u,v,w,y) = u 

D(T,T,u,v,w,y) = v 

E(T,T,u,v,w,y) = x 

F(T,T,u,v,w,y) = y. 

We have used the following notations for the characteristic functions of du j and dux 

= [ e zl dvj{l) 



(4) 

(5) 
(6) 
(7) 
(8) 
0) 

(10) 

(11) 
(12) 

(13) 

(14) 

(15) 



4>x{z) 



e zx dv x {x) . 



As a consequence, equation ^ holds if and only if these differential equations are satisfied with 
their boundary conditions. □ 

2.1.2.3 Analytical solution These equations can be solved in closed form, when parameters are 
piecewise constant. We consider first the case where all parameters are constant, which completes 
the proof of theorem [TJ 

The equations for C, D, E and F have the simple solutions 

C(t,T;u,v,w,x,y) = u 

D(t,T;u,v,w,x,y) = v 

E(t,T;u,v,w,x,y) = x 

F(t,T;u,v,w,x,y) = y. 
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The equation for B is a Riccati equation with constant coefficients: 

d t B(t, T; u, v, w, x, y) = -\<? 2 B 2 + kB + xj) (16) 

with 

i/>(u, v,x,y) = l- e>j(n) [l - £(1 - e*)] + a(l - e x ) + C(l - . 
The standard way for solving this equation is: 

1. find a particular (constant) solution b; 

2. write B = b + / , which gives a Bernoulli equation; 

3. change the variable to g = 4 to get a linear equation. 

Instead of this, we give a more pedestrian but more direct derivation of the solution, which is 
more in the spirit of the solution of the equation for A. We rewrite equation (16) as 

dB(t,T;u,v,w,x,y) _ 
-\a 2 B 2 + kB + V> 

which can be integrated directly. 

We find the simple poles decomposition after rewriting again this equation as 

dB a 2 



[B - B+)(B - B-) 2 



-dt 



where B± are in fact the constant solutions of the differential equation (16): 

K± y/K 2 + 2cr 2 V> 
B± = y . 

With non zero k and u on the imaginary axis B + and B- are distinct and the simple pole decompo- 
sition 

dB dB 2(B + -B) 

-dt 



B -B, B-B- a 



2 



= -y/K 2 + 2(T 2 %l)dt 

is integrated as 

where we simplify slightly the notations by writing B = B(t,T;u,v,w,x,y) and B? = B(T, T,u,v,w,y). 
The exponential of this equation is 



B — B— Bj* — B— 



-^^^= X {t,T ]Ul v M y) 



B — B+ Bj 1 — B^ 

where we have introduced the notation x(t, T; u, v, w, x, y). It gives the following expression for B: 

B = B - ~ XB+ =B_ + B+ ~ B - . (17) 

When T — t goes to infinity \ goes to zero and B goes to S_ . 
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We are left with the equation for A. For the moment, we set (3 = and r\ = (the last terms 
of equation Q will give simple contributions that we give at the end of this section) . This equation 
will be solved as a function of B. Dividing equation Q (without the last term) by equation ^ gives 



n\ O0 B + 1 (ct)x(B)-l) 



dA _ 

dB ~ \<j 2 B 2 -kB + (e D (j)j(C) - 1) 



The assumption we have made on the jump size distribution ensures that the right-hand side is a 
rational function. Equation ([TJ gives the characteristic function of this distribution as 

( p x (x) = (i-ezr {n+1) . 

Therefore all poles of the rational functions are known: B + and B- are simple poles and there is a 
pole of order n + 1 in 



B e 



1 



dA _ kX^B + 7 ((1 - 9B)~( n+ V - 1) 
dB ~ 



\cr\B - B+)(B - BJ) 
To write down the pole decomposition, the first step is to use 



1 



1 



\a 2 {B - B + ){B - BJ) \(J 2 (B + -B_)\B-B + B-B^ 
to separate the fraction in two parts: 



1 



F(B) F(B) 



dA _ 

d~B~\a 2 {B + -B_)\B-B+ B - B. 



(18) 



where 



F(B) = k\ 00 B- 1 + 



The first two terms give a contribution to ± B _ B ' ± 

k\ocB± - 7 



\n+l 



(B - B e ) n+1 



±ftA„o ± 



B-B± 



where the constant term will be simplified between + and — components. The other poles of ± B _ B± 



come from the term ± 



l(-B e ) n+1 
(B-B )«+i 



(B — -Bg) n+1 " ±l[ ~ Bd) 
Multiplying by B — B± and setting B = B± give 

k± 



+ E 



B-B ± (B - B e 



{B± - Be) 



n+l ' 



The decomposition is completed by computing 



1 



B - B± \(B 



n+l 



(B± - Be) 



n+l 



E 

fc=0 



(B± - B e ) n - k+1 (B - B e ) k+l 
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The pole decomposition finally reads 
F{B) 



B - B± 



± 



k\ 00 B± + 7 



(-Be) 



n+l 



(B± - B e ) n+1 
-B e ) n+1 



B — B± 



k=0 



(B± - BqY-^IB - B 6 ) k + l 



which has to be plugged into equation (Il8|). 



We can now integrate ( 18 ) analytically: 



A = At + 



1 



\& 2 {B + - B_ 
k\ 00 B± + 7 
1 



1 



(1 - 0B±) n + 1 
1 - OB 



In 



T 7" 



1 - 0B±) 



n+l 



In 



1 - 6B T 



n 1 

± 7 ^fc(l 
k=l y 



1 



(1 - 0B±) 



n-k+l 



1 



Bt — B± J 



1 



(1 - 6B) k (1 - 9B T y 



For /3 ^ or rj ^ 0, there are additional contributions to A 

-(3(l-e x )(T-t)- V (l-ey)(T-t) . 



(19) 



(20) 



Equations (17), (19) and (|20|) together with terminal conditions Bt = w and = provide 



an analytical solution for the four-dimensional characteristic function of the joint distribution of Lt, 
Nt, At, Qt and Rt- 

The formulas computed correspond to the case where all parameters are constant between the 
computation time and the maturity T for which the distributions must be computed. If parameters 
are piecewise constant, the integration has to be performed backwards, with intermediary boundary 
conditions Bt and At ■ 



2.1.3 Loss distribution and moments 

An inverse Fourier transform gives the Loss distribution from the characteristic function 



E, ( e ipLT 



e A{t,T;ip,0,0,0fi)+B(t,T;ipfi,0fi,0)\t+ipL t 

where u = ip and v = w = y = 0. 

The number of defaults can be obtained similarly by Fourier-inverting 

E ( e ipN T ^ _ e A(t,T;0,ip,0,0,0)+B{t,T;0,ip,0,0,0)\t+ipN t 



(21) 



where v = ip and u = w = y = 0. Computations involving Nt are similar to ones with Lt so in what 
follows we focus on Lt. 

We have to take into account the fact that the whole basket may default, which is modelled by 
the nonnegative stochastic variable Q. Q > means that all names have defaulted. We write 



i 



Qt=0 



-coQj 
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This allows to compute quantities like 

E t (e iplT l QT=0 ) = E t ( e ^-°° c 



and similarly 



D A(t,T;ip,0,0,-oo,0)+B(t,T;ip,0,0,-oo,0)\t+ipL t -^ 



Qt=0 



E t \e ipNT lQ T=Q ) = e A (t,T;0,ipfl-^fl)+B(t,T;0,ip,0-^,0)\ t +ipN tl ^ =o _ 

We need in addition 

U7 /n \ — 1 H7 / Tl \ — 1 „A(t,T;0,0,0,-oo,0)+S(t,T;0,0,0,-oo,0)Afti 

J^t (, I Q T >oJ = 1 — JH-t (1q t =oJ = 1 - e v ' ; >.>>>>> ; « iQ t =o • 

In fact, x = — oo means that e x is set to in formulas of theorem [T] 

The counterparty risk is handled similarly through variable Rt which is null until the counterparty 
defaults. In the following, we do not take it into account in formulas in order not to overload formulas. 
It is straightforward to extend them to this case. 

The loss mean E(Lt) or the loss variance E I ( Lt ) — El Lt ) I can be computed in two ways. 



Either the Taylor expansion of equation (21 ) in p ~ is computed, which produces all moments up 



to a given order, or differential equations are directly written and solved for the moments. 



2.1.4 Time dependency 

To be able to match market prices of index CDS at all maturities, some parameters must be time- 
dependent. There are several possible choice. The long-run intensity Aqo can be taken piecewise 
constant for example. An other possibility is a change of coordinate on the time variable, t is taken 
to be an increasing function i(r) of real time r. If time goes faster, defaults are more likely to happen 
in a given period of time, if times goes slower defaults have lower probability. For example if one 
wants to fit exactly the index CDS quotes at several maturity, the simplest solution is to take a 
piecewise affine function for /, with slope change at CDS maturities. When the slope changes by a 
factor a, this is equivalent to a scaling of parameters 



A - 


-> a 


X 


7 - 


-> a 


7 


9 - 


-> a 





n — 




n 


K — 


-> a 


K 


a — 


-> a 


a 


a — 




a 


P - 


-> a 


P 


a - 






c - 




c 


v - 


-> a 


V 



2.2 Finite basket 

The model described really corresponds to an infinite pool of issuers: there can be any number of 
defaults, even larger than the basket considered. There is no exhaustion of the issuers pool. Let us 
consider now a basket of Nm issuers. We denote by Nt the number of default in the finite basket and 
Lt the corresponding loss. 
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2.2.1 Truncation 



A first solution would be to truncate the number of defaults to the basket size Nm and the loss to 
Lm = IiNm- 

N t = mm{N t ,N M ) 
L t = mm(L t ,L M ) ■ 

The expected loss E(L t ) may be approximated to E(L t ) if the contributions from L t > Lm are 
negligible. If this approximation is not valid, a lengthier computation must be performed: a CDO 
tranche with detachment point Lm has to be priced to get the value of a index CDS. This would 
slow down our calibration process. Moreover, the probability that a given issuer defaults in a basket 
where many other issuers have already defaulted would be much greater than its default probability 
in a basket without defaults yet, for the same default intensity: the number of default per unit of 
time does not depend on the number of issuers alive in the basket. In some sense this is not a real 
problem: it introduces some default contagion. On the other hand, this default correlation is not 
controlled and does not come from the default intensity process. 



2.2.2 Subpool and uniform probability 

We therefore prefer an other solution. We consider that the finite pool is a part of the infinite pool. 
The infinite pool can be considered to model the whole economy. With uniform probability on pool 
components, the probability that a default of the infinite pool concerns the finite pool is proportional 
to its number of elements. Starting from a finite pool size of Nm, after N defaults on the finite basket 
a given default of the infinite pool has a probability proportional to Nm — N to occur in the basket. 
When the number of defaults in the basket is N = Nm there is no more default: the number of 
defaults is correctly bounded. We normalize this probability as N %~ M N ■ There is no loss of generality 
with this choice of normalization: the model would be the same if this probability is multiplied by 
some factor and the intensity process on the infinite pool is divided by the same factor. 

We forget for the moment the probability of a simultaneous default of the whole basket. We can 
compute the probability of k defaults in the credit basket as a sum over probabilities of j defaults 
in the infinite pool weighted by the probability that j defaults in the infinite pool gives k defaults in 
the finite basket. Denoting by pjk this probability we have 

oo 

¥(N T = k\Q T = 0) = Y, HNt = J\Qt = 0)p jk . (22) 

3=0 

The pkj coefficients satisfy the recursion 

k N M - k + 1 

Pj+1 * = iV^ fc + N M ^ (23) 

with initial condition 

Po,k = ho • (24) 

It is possible to find a closed form solution to this set of equations. The observation of the solution 
for k = 0, 1, 2 leads to the hypothesis that pjk, k < j, can be expanded in eigenmodes 

N M , l 



J- 1 M / 

Pjk = ^ a-kii 
1=0 ^ 



N M 
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Equation (23) is satisfied if and only if 

(I - k)a M = {N -M-k + l)o fc _i,i 



which leads to 



a-ki 
a>kl 



N M -l 
k-l 







aui-iy- 1 k>l 
k<l . 



The coefficients an are given by the initial condition ( 24 ) . The expression for pofc is triangular and can 
be computed recursively. To shorten the computation, we introduce a polynomial in some variable 
Z: 

N M 

P{Z) = Y,PvkZ k . 

k=0 

The initial condition pofc = <5fco is equivalent to P(Z) = 1. Using our solution for pofc hi terms of a k i 
this polynomial is 

N M k 



N M -l 



k-l 



k=0 1=0 

Switching sums in k and I we can sum over k and get 

N M 



a u (-l) k - l Z k . 



p(z) = J2^-z) Nm ~ 1 z 1 = i. 



1=0 



The unique solution to this equation is 



an 



N M 
I 



After some rearrangement the solution finally is 



Pjk 



N M 



(25) 



Using this solution and taking into account the disaster event modelled by Q, the probability of 



default in the finite basket (22) is 
F(N T = k) 



N M 
k 



1)*E (?)(-!)' 15 



+ 6 ktNu R(l QT>0 ) (26) 



where the expected values on the right-hand side are given by theorem [T] 



E 



Nm 



]t=0 



At 0,T,0,lnl 



( >,0 ,-oo,0)+B(0,T,0,ln( J- 



,0,-oo,o) A 



E(1q >o) = 1 — e yl ( ' T ' ' ' > _oo ' )+ B ( ' T ' ' ' '~ oo > ) A 

More generally, the conditional probability at time t is given by a slightly modified formula: 



(27) 
(28) 



>t(N T = k) = l k > Nt 



N M - N t 
k-N t 



l=N t 



k-N t 
l-N t 



-lVE t 



N M 



N T -N t 



1Qt=0 



6 k , NM E t (l QT>0 ) (29) 
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where the expected value on the right-hand side are given by theorem [T] 

~ / 7 \ N T ~N t 

E((1<3t>o) = 1 - e-*. T .».».».-~.»)+«'.^»-».a-~.»)^l <3t , . (31) 

We have finally obtained a closed formula for the default distribution as a finite sum. For small 
enough Nm or k it can be used. However for large values of Nm and k this formula is completely 
useless for numerical applications: it is a alternating sum of numbers with very large coefficients. For 
example for Nm = 125 and k = 62, which corresponds to the highest contribution to the standard 
senior tranche on CDX with a recovery rate of 40%, the highest coefficient has an order of magnitude 
of 10 54 and the final result is between and 1: the computations should be done with very high 
accuracy and cannot be done in a direct way on current processors. In the remaining parts of this 
section we therefore give other computation techniques more suitable to numerical implementation. 



,0,-oo,0 )+B(t,T,0,\n( 



,0,-oo,0 A t _ 



L Qt=0 



(30) 



2.2.3 Expected loss 

In order to compute CDS prices, the expected loss 

E(L T ) = hE(N T ) 



(32) 



can be computed from the mean jump size l\, equal to 1 minus the average recovery rate, and the 
expected number of defaults. To compute the latter quantity, we introduce the conditional mean 
value 

N M 



e 3 = E(iV T |iV T = j f Q T = 0) = Y^Pjkk 



k=0 



which does not depend on T as the probability of jumps in the basket does not depend explicitly on 
time. In a first step we condition all probabilities by Qt = 0, i.e. we do not consider a default of the 
whole infinite pool. The recursion equation (23) translates for the mean value to 



e j+ i 



EN M - k N M - c j 



k=0 



N M 



N m 



From the initial value eo = the solution is 



N M 



N m 



The unconditional expected number of defaults is obtained by averaging over the number of defaults 
j in the infinite pool and taking into account the default of the whole infinite pool: 



E(N T ) = Nm 1 — IE 



1 



N M 



Qt=0 



The expected value on the right-hand side corresponds to the characteristic function we have com- 



puted in theorem jlj with v = ln^l — ^t- 



-oo, u = w = y = 0: 



E(N 7 



N M 



3 ^0,T;0,ln(l-^),0,-oo,o)+B(o,T;0,ln(l-^;),0,--oo,o)Ao 



(33) 
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We find similarly the conditional average number of default: 



E t {N T ) =N M - {N M - N t )Et 



N M 



N T -N t 



Theorem [T] gives the explicit expression 



E t (N T ) = N M - (N M - N t )e 



A[t,T;0Ml-J^-) ,0,-oo,0 ) +B[ i,T;0,ln( 1- ^ ) ,0 ,-oo,0 ) A 



(34) 



2.2.4 Expected tranche loss 

For CDOs we have to compute the expected value of the tranche loss. We consider first the case of 
fixed recovery rate, with default size l\ for all components. From a tranche with attachment point a 
and detachment point d, the tranche loss is 

L a rf = (d-a) + (a- L T ) + - (d - L T ) + . 



The expected value of 



N M 



E[(K - L T ) + ] = ^P(iV T = k)(K - hky 



k=0 



can be rewritten using (22) as 



E[(K-L T )+]=J2 



j=0 



N M 

^PjkiK - hky 



k=0 



N T = j, Qt = 



(We suppose that the tranche does not contain the end of the basket: K < 1\Nm-) This quantity is 
the integral of the function 

N M 

Mj) = ^2MK-hk) + 

k=0 

against the probability distribution of iVy. 

Theorem [T] gives the Fourier transform of this distribution. We can invert it at each time T for 
which we need to compute the tranche loss and integrate against fx- When calibrating, these Fourier 
transforms has also to be done for all sets of model parameters. Alternatively, we can compute 
the Fourier transform of fx, which has to be done only once for each value of K and integrate it 
against the characteristic function of Nj<. This considerably reduces the number of numerical Fourier 
transform to perform. 

For each value of K needed, the Fourier transform of fx 

oo oo Nm 

kip) = eipj fx(j) = E eipj Y.p^ K - 

j=0 j=0 k=0 

is computed numerically, using FFT for example. A bound on j has to be estimated. Matrix 



elements pjk have an explicit form given in equation (25) but we have seen this expression is often 
useless numerically. We therefore compute these numbers directly using the recursion (23) from initial 



condition (24). As these numbers does not depend on the model parameters or even the attachment 



or detachment point but only on the basket size Nm, this computation has to be done only once for 
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all instruments on a basket of size Nm- The Fourier transform itself has to be done only once for 
each detachment point K, basket size Nm and recovery rate 1 — l\. 

We finally obtain E[(K - L T )+] as 



1 f+°° / ~ \ 

E[(K - L T ) + ] = —J f K (-p)E[e^l QT=0/ 

or more explicitly using theorem [T] 

EUK-Lt) + }=— / / Ar (_ p ) e A(0,T ; 0,ip,0,-oc,0)+B(0,T;0,ip,0,-oc,0)Ao _ (35) 
27T J_ 00 



For N th -to-default, we need to compute at every time step the quantity E(tN T< k) = F(Nt < k). 



rth 

It is computed similarly for k < Nm as 

"+00 



1 f +oc / 
E(liV T<fc ) = r- / 5fc(-p)E e JpiVT lQ T=0 
^ J -00 v 
+00 



1 r+00 

_ / ^ fc /_^ e A(0,T;0,ip ) 0-oo ) 0)+B(0,T;0,ip,0 ) -oo,0)Ao 



—00 



with 

00 00 k— 1 



= J] e^l iVT<fe = J] J] p iV 
i=o j=0 fe'=Q 



2.2.5 Stochastic recovery 



Stochastic recovery rate can be included in the framework. The loss process is no longer Lt = l\Nx- If 
recovery rates are independent of each other, Lt is given by the convolution of the jump distributions 
up to the Nx th jump. 

When computing tranche losses, the function fxU) has to be modified. More precisely, {K — l\k) + 
has to be replaced by the expected value of (K — Lt) + conditional to Nt = k: it is the integral of 
this function weighted by the distribution of Lt conditional to Nt = k, i.e. the convolution of 
distributions of the k first jumps. 

If all recovery rates have the same mean 1 — h, the expected loss is still given by E{Lt) = 1\E{Nt). 
If the recovery rate mean depends on the number of past jumps, it is no longer true. In many cases 
it is still possible to compute a simple formula for the expected loss. In other cases, one can always 
use the algorithm for tranche losses for a tranche covering the whole basket. 

3 Pricing 
3.1 Index CDS 

The price of a CDS is given by the difference between a protection leg which pays the amount of 
defaults on the basket^] 

CL CD s= [ cLL t ZC(0,t) 
Jo 

2 We write payoffs with integrals to simplify notations; they are in fact discrete sums. 
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and a premium leg which pays a fixed premium s 



FLcds = s [ dr(N M - N T )ZC(0, r) . 
Jo 

T is the maturity of the contract, time is the pricing date and ZC(0, T) is the zero coupon discount 
factor. The price at time is the risk-neutral expected value 

HCLcds - FLcds) = C dE(L T )ZC(0, r) - s F dr(N M - E(N T ))ZC(0, r) . (37) 

Jo Jo 



The expected values which appear in the right-hand side can be computed using equation (|34j). With 
a recovery rate of fixed mean size 1 — 1%, the expected value of L T is E(L T ) = liK(N T ). For more 



complex cases, see section 2.2.5 



3.2 CDO tranche 

For a CDO tranche of attachment point a and detachment point d, the tranche loss is 



L 



a,d 



if L T < a 

L T — a if a < L T < d 
(d — a) if L T > d . 



The credit and premium legs are 

CL 
FL 



- a,d 
CDO 

a,d 
CDO 



fT 

/ dL^' d ZC{0, t) 
Jo 

sj dr ((d- a) - ZC(0,t) . 



Introducing 



M, 



a,d 



(d-a)- L a T 



a,d 



The price is the risk-neutral expected value 

r-T 



E(CL^ d DO - FL% d DO ) = - dE^M^ ZC(0,r) - s J dr E ^M" ,d ^ ZC(0, r) 



(38) 



We decompose the remaining notional on the tranche as 



M, 



a,d 



(d - L T ) + - (a - L 7 



Its expected value is 



E (m*A =E[(d- L T )+] - E [(a - L T ) 



which can be computed using equation (35). This formula involves an integral in Fourier space. In 
fact only one integral in Fourier space is needed when computing the price of a CDO tranche: by 
linearity, all quantities in the Fourier space can be summed before doing a final integral. 
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3.3 N th -to-default 

For a fc th -to-default let us introduce 



It = ^-N T <k 



The credit and premium legs are 

CL k = -(1-R) f dI*ZC(0,r) 
Jo 

FL k = s J dr(l-/ r fe )zC(0,r) . 



'0 

The price is the risk-neutral expected value 



E(CLk — FLk) = ~h J tffi(/*) ZC(0,r) - a J dr (l -E(i*)) ZC(0,r) . (39) 
The probability 



E 



(/ r fc ) = P(iV T < fc) 



is computed using equation ( 36 ) . 



3.4 Options on Index CDS 

As the model we consider is a dynamic model, options on index CDS (swaptions) can also be priced. 
The owner of such an option get the right but not the obligation at some maturity to enter into an 
index CDS at a rate K fixed in advance. A Payer Swaption gives the right to pay the premium leg 
and get protection on the basket. A Receiver Swaption gives the right to receive the premium leg 
and sell protection on the basket. 

The price of the swap at exercise date t is a stochastic process which depends on At, Lt and/or 



N t . Expected values in formula (37) have to be evaluated at date t: 

MCLcds ~ FLcds) = [ T dE t (L T )ZC(t, r) - K C dr [N M - E t (N T )) ZC{t, r) . 

Jt Jt 

A Front End Protection clause can be included in a Payer Swaption: in this case if there are some 
defaults before the option expiry the corresponding amount is paid at the exercise. A Receiver 
Swaption would not be exercised if defaults occur prior to expiry. For a Payer swaption with Front- 
End Protection the price at exercise date t is 

[MCLcds -FL CDS )+L t } + 
whereas for a Receiver swaption it is 

IMFLcds ~ CLcds) ~L t } + . 
The price of the Payer swaption is the risk-neutral expected value 

SP = Eo{[E t (CL C DS ~ FLcds) + L t } + ) ZC(0, t) 
and similarly for the Receiver 

SR = E Q ([E t {FL CD s ~ CLcds) - L t ) + ) ZC{0,t) . 
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In all cases the price is conditional to the values of At, Lt and/or Nf. The information available 
at exercise date is supposed to be fully contained in these variables. The probability density of these 
variables is obtained in the following way. The joint probability distribution of Nt and At (conditional 
to Qt = 0) is computed by Fourier inversion from the characteristic function 



En 



jvNt+wXt 



computed using closed formulas ([2]), (17) and (19). This joint distribution is combined with the 
probabilities of Nt and Lt conditional to Nt to get the joint probability of Xt, N and Lt. With fixed 
default size h, the joint density is 



dP(A t , N t = k, L t ) = dP(A t , N t = J, Qt = 0)p jk 6(L t - hk)AL t 

3=0 

+ 6 k>NM dF(X t , Q t > 0)5{L t - hk)dL t 



The conditional price is then integrated against this probability density. In the most general case, 
a 2-dimensional Fourier transform and a 3-dimensional integral should be performed. In fact Lt and 
Nt carry almost the same information and Lt can be removed from the integration variables and 
supposed to be equal to hN. If the recovery rate is deterministic this is exact, in other cases it is 
a slight approximation. In any case, we have to compute a 2-dimensional Fourier transform and a 
2-dimensional integral. Making few approximations, the numerical computations can be reduced to 
few simple integrals, which can be useful in particular for calibration purpose. This will be detailed 
in section 13.7.41 



3.5 Options on CDO tranches and exotic products 

Options on CDO tranches can be priced exactly: the conditional price at the exercise date t is 
integrated against the probability density of At and L t . Other exotic payoffs can also be priced as 
the model allows to compute the expected value at any date of any future loss. 



3.6 Counterparty risk 

The default risk of a counterparty can be handled within our framework through the auxiliary process 



R introduced in section 2.1.1 R is set to at the pricing date. It jumps of one unit when the 
counterparty defaults. 

If a quantity will not be paid if the counterparty has defaulted, the characteristic functions of 
theorem [T] are computed with the last parameter set to y = — oo (or equivalently e v = in all 
formulas). In this limit, the contributions to the expected values are set to zero when Q is positive. 



3.7 Large pool approximation 
3.7.1 Limit 

The distribution of Nt conditional to Nt and Qt = that we have denoted by pjk = F(Nt = k\NT = 
j, Qt = 0) has significant values around its mean value 

ej = E(N T \N T = j, Qt = 0) = N m 
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The variance can also be computed: it is 

71 ■ i — t 

"3 



Vn 



1 e* = N M 



k=0 



We consider a fraction p of the portfolio, i.e. an average number of defaults ej ~ pNm- It corresponds 
to 

ln(l — p) 
3 ~ ln(l - 1/N M ) " 

At fixed portfolio percentage p the variance asymptotic when Nm goes to infinity is 

Vj -p)N M ■ 

In proportion of the basket size Nm, the standard deviation is therefore 



N M 



lp(l-p) 
N M 



When Nm goes to infinity, the distribution becomes peaked around its mean value. Our large pool 
approximation consists in neglecting the dispersion of Nt conditional to Nt- we take 



N T ^ N M 



1-1- 



1 x N T 

1Qt=0 



N M J 



(40) 



Similarly for the loss process we take 



1 - e~^ Lr l 



Qt=o 



(41) 



with Lm — 1\Nm and \i such that 



K-A*) = / « 



-,ul 



N M 



The value of the constants are chosen so that the mean of Lt conditional to Nt is 1±Nt- (By definition 
the mean default size l\ is l\ = f jduj(j).) 

This is an acceptable approximation if the basket size Nm is large enough. Indeed for low losses, 
L ~ L and jump sizes are not reduced. With many defaults, the recovery rate distribution smooths 
the distribution so that the convolution implied by this approximation does not modify it too much. 
Moreover senior CDO tranches are usually wide which means that the short scale information of 
the distribution that we lose is not very important. The only case where it would be important to 
perform an exact computation would be the pricing of a N th -to-default where the combinatorics has 
to be performed or at least approximated. As far as only the integral of the loss or default number 
on some range is considered, this approximation is reasonable. 



3.7.2 Expected loss 

In order to compute CDS prices, the expected loss 



Lm 



l-E(e-^lQ T=0 ) 
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can be computed analytically: it corresponds to the characteristic function we have computed with 
u = —fj,, x = — oo, v = w = y = 0: 



E(L T ) = L M 



1 _ e A°> T ;-M.O,0,-oo,0)+-B(0,T;-/Lt,0,0,-oo,0)Ao 



As we said above, constants Lm and fi are such that it is equal to 



E(L T ) = hE(N T ) = hN M 1 — E 



iV r 



By construction, this is in fact the exact expression of section 2.2.3 



3.7.3 Expected tranche loss 



For CDOs, the pricing of section 3.2 using exact formulas is fast enough. However we show here how 
it can be done in our large pool approximation. It is similar to the pricing of equity options with 
Levy models. We have to compute the expected value of the tranche loss. It can be obtained as the 
difference of two terms of the form 



E[(K-L T )+] =L m E 



with 



K 



1 / K 
-In 1- — 
fi V L M 



(42) 



To compute this quantity we introduce a second distribution dP(Lx, Qt) f° r the joint law of Lt 
and Qt with an explicit expression for the density function, such that it has the same mean value of 
e~^ LT \Q T= Q as the real distribution &¥{Lt,Qt) coming from the double stochastic process, i.e. the 
same mean value of Lt. 

The difference between both distributions for this quantity is computed by inverse Fourier trans- 
form from 

"+°° ~ ~ / ^\ 

dKe uK (E — E) [(K — L t ) + ] 



" k 1 (dP(L T , Q T ) - dP(L T , Q T ) ) L M {e~^ T - e^ K ) 1 Qt=0 



"+oo _ rK 

dKe 

-oo J — oo 

L M t + (dP(Z T , Qt) - dP(L T , Qt 
L M t +C ° (dP(L T , Qt) - dP(L T , Qt 



+oo 



dKe uK > 
1 1 



-iiK 



1t=0 



[iL 



M 



u — n u 



{u ~» )Lt 1q t =o 



u(u — fl) 

where we have used the fact that by construction of d£^ 

"+OO 



+oo 



dP(L T ,Q T )-dP(L T ,Q T ) 1 = 



dP(L T ,Q T )-dP(L T ,g T )) e -^lQ T=0 = 0. 
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Inverting this Fourier transform we get 



E[(K - L T ) + ] =E[(K- Lr)+] + 



1 

2^ 



+00 



dpe 



-ipK _ 



ip (ip — n) 



E 



where K is given by equation (42). Expected values E under dP(L*r, Qt) are computed in an explicit 
way. E^e(* p-/ ^ T l(2 T= o^ is given by theorem jlj 



3 40,T;ip-At,0,0,-oo,0)+S(0,T;ip-/x,0,0,-oo,0)Ao 

There are many choices for the auxiliary distribution dP(L-r, Qt)- We give two examples: 
• The simplest choice is 

^ ~ /- 1 r / ~ \i\ 

$Qt,o ■ 



dF(L T ,Q T ) = 5 1 L T + -ln\E(e-» L *t QT=0 

It corresponds to the computation of the "time value" of the tranche loss by Fourier transform, 
which is then added to the "intrinsic value" 



E[(K-L T ) + ] =L M (E 



-fiK 



• Numerically it may be a better idea to take the product of a distribution on Qt such that 

F(Q T = 0) = P(Q T = 0) 
and a Poisson distribution on with parameter 



In 



At 



E( e-^Tl 



Qt=0; 



E(1Qt=o) 



1 - e"^ 1 



For this distribution the quantity E[(if — Lt) + ] = LmE 
computed exactly from the Poisson distribution as a finite sum 



e -^L T _ e -^K\ 1( 



can be 



Numerically, one can perform the Fourier inversion either using the FFT algorithm or using a 
one-dimensional quadrature. If many detachment points K are needed the FFT is the natural choice. 
In the usual case where only one or a few tranches are priced (for example the five tranches of iTraxx) 
performing the one-dimensional integrals directly will be faster. 



3.7.4 Fast pricing of index swaptions 

Proposition 2. The price of a Receiver Swaption of strike K, maturity t and tenor T — t with Front 
End Protection can be computed (approximately) as a one- dimensional integral in Fourier space. 

Let define the function 



h(r) 



5{t-T) + ^ + r(r) 



ZC(t,r) 
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where r is the instantaneous interest rate and 

rT 



L^(A t ) = ilnQ" dr/l ( T ) e ^;-M,0,0,-oo,0)+B(i,r;- M AO,-oo,0)A^ _ 

// L±(0) is negative, the Receiver Swaption is never exercised and its price is zero. Otherwise, let 
A > be a typical valuPjfor X t such that L*(A) > 0. Introducing 



(3 A ~- 



- In (J dr ^( T ) e 4*:-r;-M.O,0 -oo,0)+B(i,r;- At ,0,0 ,-oo,0)A^ 
1 f T dr/i(r) e ^ ; -^ ' '- oo ' ) + ^ T; -^ ' '- oo '°) A B(t, r; -//, 0, 0, -oo, 0) 



T 

^(t,r;-/i,0,0,-oo,0)+B(i,T;-yU,0,0,-oo,0)A 



dr/i(r)e 

i/ie price o/ the Receiver Swaption can be approximated by 
SR = L M 1 - F° dp Sln(p(aA f ^ AA)) [ I* dr MT ) e ^ ; -M,o,o,-~>,o) 



J_oo P 

]g^ e -«p£t+ip/3AAt-/iZi+B(t,r;-/i,0,0,-oo,0)A i ^ _ jg ^ g -ipL t +i^ A At lg t _ J ZC(Q,t) 

with 

]£^ e ~ ipLt+ipPA^t— l*Lt+B(t,T;— /i,0,0,— oo,0)At^ _ 

e 40,t;-ip-^,0,ip/3 A +B(t,r;-M,0,0,-oo,0),0]+B[0,t;-ip-At,0,ip+B(i,r;-/i,0,0,-oo,0)/3A,0]Ao 

and 

jg/'g— ip£t+ip/? A At -i _ \ _ e A(0,t;-jp,0,jp/3A,-oo,0)+-B(0,i;-ip,0,ip/3A,-oo,0)Ao 



T/ie price o/ a Payer Swaption is obtained by Call-Put parity: 



SP = L M 

with 



1- / dr/i(r)E(e %% =0 



ZC(0,t) - SR 



Lr 



Ele^^ln =n I = e A0.T;-M,0,Q,-DO,Q)+fl(0,T;-^Q J 0,-oo 1 0)A 



?r=0 

As Call-Put parity holds for Receiver and Payer Swaptions, we consider here only a Receiver 
swaption. The price of a Payer Swaption is given in term of the Receiver price as 

SP = E (CL CDS - FLcds + L t ) ZC(0, t) - SR . 

Proof. First, replacing N by L/li, which is an approximation if the recovery rate is stochastic, the 
option holds on the quantity 

E t (FL C DS ~ CLcds) -L t = ^- j\r [L M - E t {L T )] ZC{t, r) - dE t {L T )ZC{t, r) - L t . 
Integrating by parts, this can be rephrased as 

MFLcds ~ CLcds) ~ L t = [ dr h ( T ) i L M - MLr)] - L M (43) 

3 Which may be as simple as A = 0. 
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where 



h(r) 



5{T-T) + ^ + r{T) 



ZC{t,r) . 



r(r) is the instantaneous interest rate: 

ZC(t,r) = e J7<M*) 
Using L T ~ Lm (1 — e _MZ/T lQ r= o), equation (43) gives 



[MFLcds ~ CLcds) ~ L t } + = L 



M 



dTh(r)E t (e 



1 



(44) 



The expected value at date t is conditional to At, Lf and Qt- According to theorem [lj it is 

e A(t,T;-^,0,0-oo,0)+B{t,T;-fifl,0,-oo,0)Xt-fJ.Lt-^^ _ q ^g-j 



If Qt = and for a given value of At, the quantity inside the bracket in equation (44) decreases with 



Lt. There is a strike value L*(At) for which this quantity becomes negative and the payoff vanishes: 



,A(t,T;-n,0,0-oo,0)+B(t,T;-fj,,0,0-oo,0)\t 



(46) 



If the interest rate is not too negative, the integral is always positive. Depending of the strike K and 
the value of At it may give a negative L*(At); in this case all values of Lf will contribute to the final 
payoff. 

In particular if L*(0) < 0, and also if Qt > 0, the Receiver swaption is never exercised and its 
price is zero: 

SR = . (47) 

In the following we consider the nontrivial case where L±(0) > 0. 

We make an approximation at this point: L*(At) is approximated by an affine function in Xt- 
Starting from a typical value for At that we denote by A, equation (46) can be rewritten as 

fT 



L*(At) = iln(7 drh(r)e 



A(t,T;-fj,fi,0,-oo,0)+B(t,T;~ti,0,0-oo,0)A 



+ ^-lnf£ dTg A (r)e 



B(t,T;-^,0,0,-oo,0)(A t -A) 



(48) 



where gA( T ) is a measure on time given by 

^ T ^ e A(t,T;-fMflfi-OCfi) + B(t,T]-^fifl-OCfl)A 



9(r) 

We finally linearize L*(Af) as 



L*(A t ) =a A + /3a (At - A) + O ((A* - A) 2 ) 



where 



OA 



Pa 



A(t,T;-n,Ofi-oo,0)+B(t,T;-Li,0,0-oofi)A 



-ln( [ dr/i(r)e 

w Wt 

- / drg A (r)5(t,T;-/i,0, 0,-oo,0) . 
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To support this approximation, we can notice first that for short tenors, the Dirac function dominates 
in h(r), which means that the measure g\(T) reduces to a Dirac delta 5(t — T). In this limit case, 
equation (48) shows that the approximation is exact. At the opposite limit, for tenors which are long 
compared to the characteristic time B(t,r,- ■ •) is almost everywhere equal to its limit value 
— -£?_(• • • )• The exponential can be factorized outside the integral, and the approximation becomes 
exact in this limit. Between these two cases (short and large tenors), the approximation can be 
controlled. The first term in the expansion in At — A which is neglected is the quadratic term. Its 
coefficient is the variance of B(t,r, • • • ) under the measure g\{r). Note that errors on L*(At) have 
no importance where the joint density of At and Lt is negligible. This is especially true for negative 
L*(At) or negative At where the density is zero. 

The price of the swaption is given by the risk-neutral expected value: the conditional payoff must 
be integrated over At, Lt and Qt using the joint distribution of these two stochastic variables. In 
order to do this, we rewrite the payoff (44) using the boundary L*(At) as 

[MFLcds ~ CLcds) ~ L t ) + = L M dr h(r)E t {e^ 1 Qt=0 ) l Z(<MAt) - ^Q t =oh t <L^x t )) ■ 

(49) 

The expected value at time of this expression will be computed by Fourier transform in Lt. Both 
terms vanish at +oo but are constant or diverge in — oo. As Lt has zero probability of being negative, 
a cut-off can safely be set at a negative value. Therefore 

1 L t <U(X t ) 

can be replaced in both terms by 

for any M > L±(0) without changing their expected value. 
The Fourier transform in Lt of this function is 

+oo ^ _ -i —ipM 

dL<e ipLt 1l ~ —- p *pMM 

OM t e L i ,(Xt)—M<Lt<L i ,(Xt) e 

-oo "¥ 

The inverse Fourier transform is the original function 

i r+oo _ -i —ipM 

t T . M<7 <T n , = — / dpe~ ipLt e^MM ( 50) 

U{Xt)-M<L t <L*{x t ) 2vr F ip K ' 

and gives the expected value today: 

Using the linear approximation L+(\ t ) ~ aj\ + f3\(\t — A), this expected value can be computed 
as a one-dimensional integral 

/ \ 1 r+oo i _ —ipM , 

E (%=olz t <^(A t) ) = ^ dp— ^e^-^ A )E(e-^^^l Qt= 
where E^e~* p ^ t+ * P ^ AA ' l<g t= o^ is computed using theorem jlj 
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As we are in the case -£/*(0) > 0, this formula can be simplified if M is set to 

M = 2{a A - f3 A A) ~ 2L*(0) > L*(0) . 

With this choice for M, the expected value is 

i p+oo 



dp 



sin(p(a A -/? A A)) E f_ ipZt+ip p AXt ^ 



This gives the expected value of the second term in formula ( 49 ) . The expected value of the first 
term 



T 



drh(r)M t (e-^t Qr=0 )l~ Lt<L , Xt 



.(At) 



is computed similarly. First, remind that the conditional expected value is given by equation (45): 



o = e 



y4(t,T;-/i,0,0,-oo,0)+B(t,T;-/i,0,Q,-oo,Q)At-/xL t j ^ 



lj t<i ( At ) is replaced by 1 L (\ t )-M<Z t <L (A t ) wmcn can be rewritten using equation (50). This gives 



E 



J 1 L4X t )-M<L t <L- t (\ t ) 



e A(t,T;-ft,0,0,-oo,0)+B(t,r,-v,0,0,-oo,0)\t-iiLt j q / dpe~ ipL ' 

27T ' ' 



+oo 



~ i _ e -ipM 



,ipL*(A t ) 



The expected value at time of this expression gives 
E[E t (e-^l Qr=0 ) lz t <^(A t )' 

1 />+oo -i _ —ipM . _ _ 
e ^(t,r;-M,0 ) Q ) -oo ) Q) J_ / d 1 e £ / e -ipL t +ipL*(A t )- A tL t +B(t,r;-At 1 0,0 ,-oo,0)A t j ( 

Using again the linear approximation L+(\ t ) ~ «a + /3a(A( — A), this expression can also be 
evaluated as a one- dimensional integral 



E 



E t f e -"M Qr=0 ) l Zt<MAt) 



a A(t,T;-n,0,0,-oofi) 



1 



dp ■ 



1 _ e -ipM 
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_ e ip(aA-PAA)TQ^ e -ipLt+ipPA^t-fiLt+B(t,T-,-fi,0,0,-co,0)\t |^ _ q 

where E(V ipit+ip ^ AA '^ Lt+ - B ( t ' T; -^ ' -°°<°) A ' l Qt=0 ) is given by theorem [l| 



]g| e —ipLt+ip[3A\t—tiIjt+B(t,T;—n, 0,0, -oo, 0)A t -n _^ 



3 40,t;-ip-At,0,ip/3A+-B(t,r;-^,0,0,-oo,0),0]+B[0,t;-ip-At,0,ip+B(t,T;-/i,0,0,-oo,0)/3A,0]Ao 



With the choice M = 2(«a — /?aA) the integral is 



E 



: "' ^ o) l Zt<L4At) 



= 4*,r;-^,0,0,-oo,0) 



7T 



dp 



sin(p(a A - /3aA)) t 
P 



The summation over coupon dates finishes the computation of the Swaption price. 



□ 
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4 Calibration 



4.1 Algorithm 

The model can be calibrated on many instruments, depending on what it is used for. As they are 
the most liquid products, the market spread of index CDS should be matched. Other relatively 
liquid instruments, at least on traded indices, are standard CDO tranches. One may therefore want 
to match CDO market quotes for all available maturities. As the model for L can be the sum of 
many basic building blocks with their own piecewise constant parameters sets, in principle many 
(non-arbitrageable) sets of price can be reached. However, introducing too many parameters make 
the calibration more difficult and more unstable. With too many parameters, the model can even be 
under-specified. So we prefer to stick with one single building block. Moreover, time-dependency is 
introduced in the minimal way of section [2 .1.4 through a piecewise affine change of time. We calibrate 



the model to match as closely as possible market quotes for CDO tranches (either one maturity or 
all maturities) while matching precisely quotes for the index CDS. (It would be useful to calibrate 
also on swaptions: CDS or CDO tranches are not really sensitive to some dynamical properties of 
the model.) 

The calibration routine uses the following steps: 

1. Choose values for parameters Ao, Aqo, k, a, n, 6, 7, a and /?. 

2. Conditionally to this set of parameters, calibrate the slopes a% of the time-dependence function 
t = t{r) by a bootstrap method. One dichotomy or Newton algorithm for each CDS maturity 
is needed, which is fast because there is a closed formula for the expected loss. CDS prices are 



computed as explained in section 3.1 



3. Price CDO tranches (as in section 3.2) for this set of parameters and time-dependency function 



and compute the sum of squared errors, relatively to the bid-ask spreads. 
4. Loop on step 2 according to some optimization algorithm. 

We use differential evolution [SP95] algorithm to find a global minimum of the objective function. 



4.2 Numerical results 

We take as an example iTraxx Europe Series 9 Version 1 (RED Code: 2I666VAI6) data on September 
30th, 2009. Quotes for the index, 9%-12% and 12%-22% tranches are spreads in bps. 0%-3%, 3%-6% 
and 6%-9% tranches have a running spread of 500bps and the quote are upfronts in percent. 



4.2.1 Global calibration 

Table [l] and figure [|0 show market quotes and the result of a global calibration on 5Y, 7Y and 10Y 
instruments (maturing on June 20th, 2013, 2015 and 2018). We set the recovery rate to 40%. 



4 We plot the quantity ^J^j 2 > wnere s is the spread quote (tranches 4 and 5) or the running spread (5% for tranches 
1, 2 and 3), T is the maturity in years and U the upfront quote for tranches 1, 2 and 3. Forgetting about discount 
factors and with the hypothesis that on a given tranche the loss per unit of time is constant, we can write the upfront 



in term of the tranche loss L and running spread s as U ~ L - sT (1 - f ) . This gives an approximation of the tranche 

lns«- T ^ sT+U 
IOSS. 1j ~ 1+sT/2 - 
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5Y 


iviarKet 
7Y 


10Y 


5Y 


iviouei 
7Y 


10Y 


Index 


1 D2 505% 


1 03 487% 


1 04 985% 


1 02 505% 


1 03 487% 

IUu.tU i / U 


1 04 985% 




165 


170 


175 








0%-3% 


36.81% (1.06) 


45.81% (1.06) 


52.50% (1.06) 


36.10% 


45.52% 


52.94% 




trnn 


ouu 


oUU 








3%-6% 


2.83% (1.06) 


8.06% (1.12) 


14.88% (1.12) 


1.49% 


8.55% 


14.93% 




500 


500 


500 








6%-9% 


-6.95% (0.88) 


-5.59% (0.88) 


-1.94% (0.88) 


-7.44% 


-5.13% 


-2.33% 




500 


500 


500 








9%-12% 


147.75 (12) 


196.13 (12) 


247.00 (11.5) 


127.20 


200.01 


259.85 


12%-22% 


58.75 (8) 


81.88 (8) 


98.88 (8) 


54.92 


80.48 


105.07 



Table 1: Market quotes for iTraxx Europe Series 9, index and standard tranches on September 30, 
2009 and model results after a global calibration. Index quotes are given as 100%— upfront, the three 
first tranches are quoted upfront in percent; the running spreads in bp are given in small figures below 
quotes. The last two tranches have quotes in spread, in bp. The numbers in parentheses are bid-ask 
spreads. 




Figure 1: Comparison of market and model values after a global calibration for iTraxx Europe Series 9 
standard tranches on September 30, 2009. The value plotted is the quantity ppgj^g , in percent, where 
s is the spread quote in percent (tranches 4 and 5) or the running spread (5% for tranches 1, 2 and 
3), T is the maturity in years and U the upfront quote in percent for tranches 1, 2 and 3. 



2N 



The parameters obtained by a global calibration on all maturities are the following: 



A 7.115 

Aoo 0.1846 

k 4.303 

a 0.9085 



n 9 

7 0.08774 

6 5.704 

a 1.377 10 _ 
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P 0.005865 
a\ 0.766 
a 2 1.268 
a 3 1.111 . 



ax, a% and 03 are the slopes of the piecewise affine function t = t(r) used to match index CDS quotes. 
Calibrating eight free parameters plus three slopes we manage to fit exactly the three index quotes 
and get a reasonable fit of the 15 CDO market quotes. Seven CDO quotes fall precisely within the 
bid-ask spread and the other are a little outside. 

4.2.2 Single maturity calibration 

In order to test the model ability to match all quotes for a given maturity, we keep five free parameters 
and freeze three quantities to the values obtained in the global calibration: ^ = 0.04289, = 
0.5195 and a = 0. We calibrate on CDO tranches of one single maturity but we still match the index 
quotes of the three maturities using the piecewise affine time change. The results are shown in table 
[2] and figure [2} All quotes are perfectly matched within the bid-ask spread. The model parameters 
for each maturity are the following: 



5Y 7Y 10Y 

A 1.013 1.815 6.379 
Aoo 0.01748 0.03061 0.3463 



K 0.4076 0.7135 8.075 
a 0.06084 0.1065 1.205 



4 17 34 



7 0.1049 0.09124 0.07845 

9 1.622 0.6226 3.131 

a 

p 0.004045 0.005947 0.006013 



01 0.996 0.851 0.834 
a 2 1.104 1.193 1.215 
a 3 0.936 1.056 1.063 
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5Y 


Market 
7Y 


10Y 


5Y 


Model 
7Y 


10Y 


Index 


102.505% 


103.487% 


104.985% 


102.505% 


103.487% 


104.985% 




165 


170 


175 








0%-3% 


36.81% (1.06) 


45.81% (1.06) 


52.50% (1.06) 


36.78% 


45.80% 


52.52% 




500 


500 


500 








3%-6% 


2.83% (1.06) 


8.06% (1.12) 


14.88% (1.12) 


2.79% 


8.05% 


14.90% 




500 


500 


500 








6%-9% 


-6.95% (0.88) 


-5.59% (0.88) 


-1.94% (0.88) 


-6.98% 


-5.59% 


-1.92% 




500 


500 


500 








9%-12% 


147.75 (12) 


196.13 (12) 


247.00 (11.5) 


148.44 


196.29 


245.97 


12%-22% 


58.75 (8) 


81.88 (8) 


98.88 (8) 


57.61 


81.40 


100.31 



Table 2: Market quotes for iTraxx Europe Series 9, index and standard tranches on September 30, 2009 
and model results after one calibration for each maturity. Index quotes are given as 100%— upfront, 
the three first tranches are quoted upfront in percent; the running spreads in bp are given in small 
figures below quotes. The last two tranches have quotes in spread, in bp. The numbers in parentheses 
are bid-ask spreads. 
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Market 7Y - K - 
Market 10Y - -* - 
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Model 7Y — ■- 
Model 10Y - o- - 




Figure 2: Comparison of market and model values after one calibration for each maturity for iTraxx 
Europe Series 9 standard tranches on September 30, 2009. The value plotted is the quantity ppg^ , 
in percent, where s is the spread quote in percent (tranches 4 and 5) or the running spread (5% for 
tranches 1, 2 and 3), T is the maturity in years and U the upfront quote in percent for tranches 1, 2 
and 3. 
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5 Conclusion 



We have presented a model for credit index derivatives with the following properties: dynamic, 
with volatility and jumps on the loss intensity, allowing closed or semi-closed pricing of CDS, CDO 
tranches, N th -to-default and options. It is also possible to price options on tranches or other exotic 
products since forward conditional distributions can be computed. 

In addition to default risk, our model handles volatility of the index spread and jump events, 
which model possible crises. Beyond pricing, this can be useful for risk management. Counterparty 
risk can also be included. 

It is possible to match exactly index prices while matching reasonably CDO tranches quotes. 
If only one maturity is considered for CDOs, we obtained a perfect calibration within the bid- ask 
spread. It would be interesting to calibrate jointly on option prices to fix the remaining parameters 
(mean reversion and volatility in particular) and get robust sets of parameters. 

It would be interesting to study our model in the spirit of |C K08| to see how it performs with 
respect to hedging of default risk and spread risk: it contains several properties that the authors find 
desirable, in particular the presence of jumps in the spread, not caused by the default of one of the 
constituents. 



31 



References 



[AH07] Matthias Arnsdorf and Igor Halperin. BSLP: Markovian bivariate spread-loss model for 



portfolio credit derivatives. 2007. http://arxiv.org/abs/0901.3398 



[BPT07] Damiano Brigo, Andrea Pallavicini, and Roberto Torresetti. Calibration of CDO tranches 



with the dynamical generalized-Poisson loss model. Risk Magazine, 5:70-75, 2007. |http : 
|//ssrn. com/abstract=900549| 

[CK08] Rama Cont and Yu Hang Kan. Dynamic hedging of portfolio credit derivatives. 2008. 



http : //ssrn . com/abstract=1349847 



[DG01] Darrell Duffle and Nicolae Garleanu. Risk and valuation of collateralized debt obligations. 
Financial Analysts Journal, 57(l):41-59, 2001. |http : //ssrn . com/ abstr act=265 223| 

[DPS00] Darrell Duffle, Jun Pan, and Kenneth J. Singleton. Transform analysis and asset pric- 



ing for affine jump-diffusions. Econometrica, 68(6):1343-1376, 2000. http://ssrn.com/ 
labstract=157733l 

[EGG09] Eymen Errais, Kay Giesecke, and Lisa R. Goldberg. Affine point processes and portfolio 
credit risk. 2009. |http : //ssrn . com/abstract=908045| 

[GGD09] Kay Giesecke, Lisa R. Goldberg, and Xiaowei Ding. A top-down approach to multi-name 
credit. 2009. http: //ssrn. com/ abstract=678966. 

[LiOO] David X. Li. On default correlation: a copula function approach. Journal of Fixed Income, 
9(4):43-54, 2000. 

[Sch06] Philipp J. Schonbucher. Portfolio losses and the term structure of loss transition rates: 
A new methodology for the pricing of portfolio credit derivatives. 2006. http: //www. 
schonbucher . de/papers/ cdo_loss_transition_rates . pdf , 



[SP95] Rainer Storn and Kenneth Price. Differential evolution — a simple and efficient adaptive 
scheme for global optimization over continuous spaces. Technical Report TR-95-012, ICSI, 
March 1995. 



32 



